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Abstract. Let a pure state I?/') be chosen randomly in an A'^Af-dimensional Hilbert space, and consider the 
reduced density matrix pA of an A^'-dimensional subsystem. The bipartite entanglement properties of \ip) 
are encoded in the spectrum of pA- By means of a saddle point method and using a "Coulomb gas" model 
for the eigenvalues, we obtain the typical spectrum of reduced density matrices. We consider the cases of 
an unbiased ensemble of pure states and of a fixed value of the purity. We finally obtain the eigenvalue 
distribution by using a statistical mechanics approach based on the introduction of a partition function. 
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PACS. 03.67.Mn Entanglement characterizations - 02.50.Sk Multivariate analysis - 05.70.Fh Phase tran- 
sitions 

^ 1 Introduction 

Li the last years, many efforts have been directed towards the study of random quantum correlations [1112] • This work 
is a short overview of some important results on the distribution of entanglement among two subparts of a large 
I '~' I quantum system. This has been a subject of interest among physicists and mathematicians for a long time, and many 
interesting results have been achieved, for instance in the context of quantum maps [3l|4]. Our presentation has the 
T-H style of a pedagogical review, in the sense that many important results will often be rederived by simpler methods. 
^ The mathematical techniques will be easy to follow, and will often be guided by physical insight. 

Our aim is to find the maximum of the joint distribution of the eigenvalues of the reduced density matrix. We 
shall therefore focus on the most probable spectrum of the reduced density matrices. In order to solve the saddle 
point equations, where the gradient vanishes, we shall invoke the physical interpretation of the so-called "Coulomb 
gas" [5], displaying the profound link between our problem and constrained 2I?-electrostatic models. It turns out 
that the most probable eigenvalues of the density matrix coincide with the equilibrium positions of movable charges 
on a line when the interaction forces arise from a logarithmic potential. It is very fascinating to discover that this 
j equilibrium configuration (the A^-tuple of eigenvalues that maximize probability) is reached at the zeros of a class of 
. . orthogonal polynomials. Therefore, the spectrum of a typical quantum state is completely determined by the zeros of 
a certain polynomial. Even for high-degree polynomials, we shall obtain some useful analytic results (in a compact 
and manageable form) for some entanglement quantifiers, such as purity. 

All these results agree with known results about the statistical averages of typical states obtained with other 
methods. We shall also recover many thermodynamic limiting results, generally obtained in literature with methods 
based on statistical mechanics or random matrix theory [^[TIISIIMTU] (for a review of random matrix theory, see Ref. 
We will try to avoid unnecessary mathematical details or unclear physical hypotheses. 

1.1 Notation: setting up the problem 

Let us consider a bipartite quantum system S ~ A -\- B whose associated Hilbert space is a tensor product 

Hs^'Ha<E)Hb, with dimHA^ N <dimnB ^ M. (1) 

If the global system 5* is described by a pure state, that is a unit vector |?/') G Tis, subsystem A is described by the 
reduced density matrix, obtained by tracing out subsystem B 

PA = tlB \^} (V-l , (2) 
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which is known to be a Hermitian, positive, unit-trace N x N matrix. In terms of its eigenvalues, the spectrum of pA 
is a probabiUty vector A = (Ai, . . . , A^r), formed by N nonnegative numbers Xj > that sum up to 1, i.e. ^ A^ = 1. 
The pure quantum states of a bipartite system "Ha ®'Hb are said to be separable if they admit the product form 
IV') = 10)^ ® Ix)b for some \4')aAx)b belonging to "Ha and 'Hs, respectively. If the state does not admit such a 
factorization, it is said to be entangled. 

Bipartite pure states represent one of the few cases for which the problem of revealing quantum correlations admits 
an exhaustive answer. From the Schmidt decomposition, one has to look at the spectrum of the reduced density matrix 
Pa', if its spectrum is (a permutation of) A = (1,0,..., 0), that is p\ = pA is a rank one projection, then the state 
is separable. Otherwise, if rankp^i > 1, the state is entangled. 

A first amount of information about the separability of a state with respect to a given bipartition is provided by 
the so-called Schmidt number |12j . defined as the rank of the reduced density matrices of its subparts. In order to 
better quantify entanglement, one can recall the statistical interpretation of the reduced density matrix: the more 
entangled the global pure state, the more mixed the reduced state. The eigenvalues of the reduced density matrix form 
a probability vector and the degree of mixedness of a state is related to how "nearly equal" the eigenvalues are. To 
this end, it is useful to have in mind the basic ideas of majorization theory |13j . We recall that the probability vector 
X is said to be majorized by the probability vector y, in symbols a; ^ y, if 

fc fc 

max^ a;^(j) < max^ yg-Qj, l<k<N, (3) 

where cr is a permutation. Thus, the eigenvalues of density matrices satisfy 

(1/7V, . . . , 1/iV) -< (Ai,...,AAr) -< (1,0,. ..,0), (4) 

in agreement with the degree of mixedness, and the corresponding degree of entanglement. 

A manageable measure of entanglement is the local purity: Given a bipartite pure state |V') with reduced density 
matrices pa and ps, one defines 

N 

TTAB = tr^ Pa = trs p| = ^ A^, 1/iV < ttab < 1 ■ (5) 

i=i 

The local purity has the Schur- convexity property, i.e., it preserves majorization order: 

TTAB (A) < TTABifJ'), if A -< /X. (6) 

As such, it is a good entanglement measure. In particular ttab = 1 and ttab = ^/N iff is, respectively, separable 
and maximally entangled with respect to the given bipartition. 

Let us start to discuss the typical properties of random pure states, that is unit vectors "drawn at random" from the 
Hilbert space. The rationale behind the sampling criterion is to introduce no bias in our judgments, namely information 
that we do not have. The idea of randomly taking a pure quantum state is then equivalent to assuming minimal a 
priori knowledge about the system. Identifying minimal knowledge with maximal symmetry, it is natural to require 
that the statistical ensemble be invariant under the full group of unitary transformations. Thus the sampling criterion 
corresponds to a unique "natural" measure on states, induced by the Haar probability measure dpniU) on the unitary 
group U{L), where L = MN = dimHs- In other words, a random pure state, defined by the action — U l^po) of a 
random unitary matrix U on a given reference state \ipo), can be represented, in an arbitrary basis, as a given column 
of the random unitary matrix U. Such an ensemble may be identified as the "most random" ensemble of possible 
states of the system. By tracing over subsystem B, this measure translates into a measure over the space of Hermitian, 
positive matrices of unit trace |14j . By standard methods one can show that the measure dp{p) of density matrices 
p G T>{N) factorizes into a product measure, with respect to the diagonalization p = J7^diag(A)[/, where U G U{N) 
and diag(A) is the diagonal matrix with diagonal entries A, 

dp{p) = diyiX) X dpniU) , (7) 

where the first factor defines a measure on the (N — l)-dimensional simplex of the probability vectors of eigenvalues 

Zi^^_l={AeM^|Afe>0,^Afe = l}. (8) 

fc 

The second factor pH on the space of unitary matrices U{N) is responsible for the choice of the eigenvectors of p. A 
unitarily invariant measure over the pure states of a composite system induces a measure over the eigenvectors of the 
reduced density matrix pA which is still rotationally invariant, i.e. a Haar measure: Think for example of the measure 
on an equator induced by a uniform measure on a hypersphere. Observe that the space of Hermitian matrices is not 
compact, while the space of states 'D{N) is the product of two compact spaces, a simplex and a sphere. 
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2 Joint distribution of the eigenvalues 

Since the information about the separabihty of a bipartite pure state \'ip) is completely encoded in the spectrum of its 
reduced density matrix pA, we will focus our attention on the typical properties of the eigenvalues Xoi pA- For random 

pure states sampled from the unbiased ensemble ^{C/|V'o) u£U{l) '^A'-hC^)): the eigenvalues of the reduced 

A''— dimensional density matrix pA = tr b{U l-ipo) (V'ol U^) a-i'e distributed according to the measure di/(A) = fM,A4WdX, 
with joint probability density function (pdf) [T^fTB] 

l<j<k<N 1<'<JV 

{NM ^1)1 

"^"'"^ - n,<,<.(M-,)!(iV -, + !)! 

is a normalization factor, assuring that i^(Z\jv-i) = ^ fN.M{X.)dX = 1. 

Let us summarize some known results about unbiased random states and the eigenvalues distribution (|9|. From 
the permutation invariance of the eigenvalues joint pdf and the unit trace condition, one immediately obtains: 

--J2{K)^N{X,) =^ (^«) = ^' V1<*<7V, (11) 

i 

where (•) stands for the expectation value with respect to the pdf ([9]). Calculation of the second moment needs more 
work. Lubkin [IT] calculated 



where [2] 





^^l/iV2\l/2 

MN+1 



(12) 



A remarkable fact is that, in order to perform his calculation, Lubkin did not use the above joint pdf but rather a 
geometric method requiring averages over the real 2A'^M-dimensional unit sphere (the joint distribution was discovered 
many years later [15 ). By denoting p s.t. pN = M — N, in the large N limit, with the ratio M/N = 1 + p finite and 
fixed, the width of the distribution becomes 

CTims ^ I ■ (13) 

Various aspects of the entanglement properties of random pure states have been studied in previous articles. 



From (11), (12) and permutation invariance one obtains the average value of purity 



k^AB) = ( E ) = ^ (^') = ^'-™s + ^ (A.)^ = . (14) 



For M/N finite and fixed, and large N one has 



(-ab)-^^. (15) 



A balanced [p — 0) bipartite large system has typical purity {ttab) ^ 2/iV [B]. 

Another measure of bipartite entanglement for pure states is the local von Neumann entropy Sab = trjp^i \npA\- 
Its average over random pure states is given by 

k=M+l 

as conjectured by Page in his pioneering work [TS', and proved in [T5]. See also [19] . 

As a final quantity that provides information on the degree of mixedness of pA, and thus of the entanglement 
of a pure state, we recall the elementary symmetric polynomials Sfe(A) — J2ji<---<jf, ^ji ' ' ' Ajt, with 1 < fc < A^. In 
particular, the elementary invariant sat (A) = Ai • • • Aat is nothing but the determinant of the density matrix det pA- It 
is a bounded Schur-concave function of the spectrum of the reduced density matrix. If the state is sampled according 
to the unbiased ensemble, the form of the moments of det pA is a straightforward consequence of Eq. ^ 
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Fig. 1. Contour plots of the function F2,m for a single qubit (A'^ = 2) for different dimensions of tlie environment M; the 
domain is the full square Ai, A2 G [0, 1] x [0, 1]. Below, the restrictions of F2,m to the simplex (i.e. the diagonal of the square 
from (1, 0) to (0, 1) ). When M > N, two minima appear in the interior of the simplex. As M increases, the repulsion between 
the eigenvalues becomes dominated by a "uniforming" effect of the environment. 



3 Most probable distribution 



In the previous section we have presented some important results about the expectation values of some interesting 
entanglement quantifiers. An alternative approach relies on the study of the typical properties. The idea is that, given 
a function of random states hQip)), vifhich is therefore itself a random variable, its most probable value h is "close" to 
its average {h). Then, the typical properties of the reduced state pA depend on its typical spectrum. We look for the 
most probable eigenvalues, that is the point (s) A that maximizes the pdf ^ on the simplex Z\jv_i. If one writes the 
function (|9| as a "Boltzmann factor" 

/nmW = Cn,m exp (-Fjv,Af (A)) , (18) 
then the points that maximize the probability / are points of minimum for the "energy" F: 

Fn,mW = - In ^^^i^ = -2 V In |A, - A, | - (M - A^) V In A; (19) 

on the simplex A^-i. 

The function -FV_m is smooth in the subset of the simplex A^^i defined by the inequalities |Ai — Aj| > and 
\i > if M > N. Moreover, when A^ — Aj — > for some i ^ j the energy diverges, Fj^^m +00. When M > N the 
same thing happens when Ai J, for some i. Thus, in the unbalanced case there exists a finite minimum of the energy 
in the interior of the simplex, and the point of minimum is a critical point of the energy, i.e. its gradient vanishes. 
For A S An^i one obviously gets \Xi — Aj| < 1 and Ai < 1 for any whence F^^m > 0, and the minimum value is 
nonnegative. 

On the other hand, when M — N, the energy function no longer diverges on the boundary, and the minimum energy 
can be attained on the boundary, and in general it is not a critical value. For example, for a qubit, i.e. N = M = 2, 
the minima are at the boundary points (0, 1) and (1,0), where ^2,2 = and VF2^2 7^ (see Fig. [I]). In fact, we will 
show that this is always the case: one eigenvalue, say A^r vanishes at the minimum point, and the problem is reduced 
to the minimization of F/v^Ar(Ai, . . . , \n~i, 0) — F/v_i^Ar+i(Ai, . . . , AAr_i). However, this is nothing but the energy of 
the unbalanced problem on the simplex An^2, whose minimum is a critical point. 

Summarizing, the problem of maximizing the pdf fN,M reduces to the problem of finding the critical values of the 
nonnegative energy function Fn^m, when M > N (or i^7v-i,Af+i, when M = N), in the interior of the simplex Z\jv-i 
(or An-2)- 
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4 Unbiased Pure States 

We will start from unbiased states, i.e. states sampled according to the unitarily invariant Haar measure. In particular, 
we are concerned with the Schmidt coefficients with respect to a given bipartition. These coefficients give information 
about the degree of mixedness of the reduced density matrices of the subparts of the global system. We are interested in 
the typical entanglement of a small subsystems of a large random pure state. Our approach will rely on a saddle point 
method: given the joint distribution of the eigenvalues of the reduced density matrix pA € VlN), we will search the 
most probable spectrum, that is, the density matrix (up to local unitaries U G U{N)) that maximizes the probability. 

The saddle point problem for typical states can be reduced to the problem of finding the equilibrium configurations 
of a system of identical movable charges on a line (electrostatic models) . These problems are elegantly connected with 
the theory of orthogonal polynomials, as Stieltjes first showed 20, ,21,. 

For unbiased states, one is able to fully solve the problem. The complete solution of the saddle point method will be 
provided for all possible (unbalanced) bipartitions. This result is the starting point to compute all quantities of interest. 
For some well-known quantities, such as purity or elementary symmetric invariants, we are able to give compact and 
manageable analytic expressions, for all N and M. Moreover, we will present typical entanglement properties for the 
unbalanced and balanced bipartition in the large sizes limit TV, M — > +00. 



4.1 Coulomb Gas 

To get a clearer insight of the joint pdf of the eigenvalues, Eq. (9]), one can invoke the physical picture of a "Coulomb 
gas" of N repelling electric charges on a segment [5]. Indeed, according to the discussion of section [sj the most 
probable distribution is the result of a constrained minimization problem for an energy function i^jv.A/, with suitable 
N < M, that can be handled by using the method of Lagrange multiplier. Namely, one has to find the minima of the 
[N + l)-variable function: 

i^tot (A, = -2 E log - I - - ^) E log - - E ^^)' (20) 

i<j I i 

that can be interpreted as the potential energy of a gas of N point charges at positions A^'s. The potential energy is 
given by the mutual repulsion of these charges, plus a part given by an external field 

K™tuai(a;,y) = -21n|a; - y\ , (21) 
(fcxtix) = -alnx + £,x , (22) 

where we have denoted a = M — N. The external potential is plotted in Fig. [2] In other words, the total electrostatic 
energy is 

£^tot(A;e) =EE^«i"t"-i(^»'-^J") + E'^-t(-^») • (23) 

i j^i i 

We are interested in the stationary points of this energy. By deriving £'tot(A;f) with respect to both the A,;'s and ^, 
we get + 1 saddle point equations: 



1 M-iV , 

E^^-1- 



(24) 



In the framework of the electrostatic model, the saddle point equations are nothing but static equations of balance 
of the forces (the derivatives of the potential energy, VxEtot — 0), with the additional constraint that the charges 
average position be equal to 1/N. 

Before trying to write down a solution for the saddle point equation, it is useful to give a picture of what will happen. 
The N charges interact via a logarithmic 2D-Coulomb repulsion, that is the form of Gauss law in two dimensions. 
Therefore, in absence of an external field, the charges will move as far apart as possible. However, there are two 
external forces acting on them. The first one is due to a charge a, with a — M — N > 0, at position a; = that repels 
the unit charges (the eigenvalues) via a logarithmic potential. This repulsion constrains the charges on the positive 
half- line. The second field is a constant force (the gradient of the linear potential ^x). In order to ensure the existence 
of an equilibrium configuration, the external potential must have a minimum. Then, ^ must be positive (necessary 
condition to have a convex potential). 
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Fig. 2. Left: Coulomb gas charges experience an external force — 9a;<^cxt(a;). The external potential (221 is the sum of a 
logarithmic part, due to a charge a = M — A'' at the origin, and a linear part governed by the Lagrange multiplier ^. Right: the 
equilibrium positions of the charges lie at the zeros of a generalized Laguerre polynomial. 



Notice that in the balanced case, a = M — N = there is no logarithmic repulsion from the origin. Therefore, one 
of the N charges will sit at the origin and repel the remaining — 1 charges with a potential — 21nx. This is nothing 
but the electrostatic problem generated by the energy function ^V_i^jv+ij thus proving that in the balanced case the 
minima are located on the boundary of the simplex A^-i, as discussed in section |3] 

The minimum of the external potential is located at Xc = oi/^. Since we expect that, in the typical case, all 
eigenvalues be located near the maximally mixed value, this critical point has to be close to l/N . By setting the trial 
value Xc — i/N, we can immediately guess that ^ ~ NM in the large-M limit. In fact we can do better and compute 
exactly the Lagrange multiplier by means of a nice trick. Observe that, for any continuous function /i, the equality 



y-v y-v HXj) _ ^ h{Xi) - h{Xj) 

Xi Xj 



i j^i 



Xi — A, 



holds. Consider now the first N saddle point equations (24 1 and multiply each of them by Xi to obtain 



2E 



X, 



X, - X. 



+ (7\f - iV) - = , \<i<N. 



(25) 



(26) 



By summing over z, and using (251 and the unit-trace condition, one finds 

N(N - 1) + N{M - iV) - C = , 



(27) 
(28) 



whence the sought multiplier reads 

C = iV(A/-l), 

as expected by the above argument. 

Once the Lagrange multiplier is known, we can obtain the typical purity "kab with the same trick by multiplying 
the first N saddle point equations (|24|) by Af and summing over i. Since 



A? -A? 



Ai — A, 

i<3 ■' 



2 ^(A, + A,) = 2{N - 1) 5^ A, = 2{N - 1), 



i<] 



we get 
whence 



2(7V-l) + (Af-iV)-e7rAB =0, 
iV + Af - 2 

T^AB 



(29) 



(30) 



(31) 



N{M - 1) 

As a check that the stationary points of the energy i?tot(-^;'?) are in fact minima, one can look at the Hessian 



matrix 



92 



dXidXj 



i?tot(A;e) 



(32) 
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It is easy to see that 



M -N 



l<i,j < N 



(33) 



Thus the Hessian if is a strictly diagonally dominant symmetric matrix with positive diagonal elements. Therefore, 
it is positive definite everywhere, and so every stationary point of Eq. (23 1 is a local minimum. 



4.2 The Solution 



The saddle point equations (24) can be tackled by using an ingenious method due to Stieltjes, that deals with 



the electrostatic interpretation of the zeros of some families of orthogonal polynomials. The first N saddle point 



equations ( 24 1 are "equivalent" to a single polynomial equation. Indeed, one can write 



2V^-^ (34) 

where 

g{x)=l[{x-Xk) (35) 

fc 

is the nodal polynomial whose zeros are the A^'s. The above identity is commonly known as Stieltjes's trick. It can be 
easily derived by noting that 

g'ix) = Hix - A.) Y: = - ^'^O E ^ E (36) 

k i ^ k i ^ j^i 

for X different from every Ai, and by taking the limit x — A^. 
Thus, the saddle point equations can be written as 

X^g"{K) + {a- iK)g'{K) = , l<i<N . (37) 

The above equations mean that the iV-degree polynomial 

xg"{x) + {a-£,x)g'[x) (38) 

has its N zeros at (Ai, . . . , Ajv), and then it is a multiple of g{x). Comparing the leading coefficients of the two 
polynomials we get that the proportionality constant is —N^, and we arrive at the differential equation 

xg"{x) + {a- ^x)g'{x) + ^Ng^x) = , (39) 

where we recall that g{x) is a polynomial whose simple zeros are the unknown A^'s. The polynomial solution of the 
previous ordinary differential equation is an associated Laguerre polynomial |21j 

4"^'^(e^) = Ec,(-x)^ c, = ^(^"^), a = M-N. (40) 

See Fig. [2j The value of the Lagrange multiplier ^ is fixed by the trace condition. Since the coefficients of x^~^ and 
a;^ are related by cn-i — cn X^i ^i, one should have cm-i — cn, that reads (M — 1)^^^^/ {N — 1)! = /N\, yielding 
Eq. ([28]). 

Notice finally that the Lagrange multiplier ^ is related to the trace of . Indeed by just summing the saddle point 



equations and using (p| with h = l, one easily obtain ^ = Tip;^^ {M - N)/N and then Trp;^^ = N'^{M - 1)/ {M - N). 
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Fig. 3. Most probable eigenvalues of a four-level system A'^ = 4 when dimHs = M > 4 varies. In the Coulomb gas framework, 
the four charges reach their equilibrium position in the external field v'ext ~ — a Ina; + ^a;, where a — AI — 4 and ^ = 4(M — 1), 
see Eq. (281. As M increases the state tends toward the maximally mixed one, i.e. A = (1/4, 1/4, 1/4, 1/4). 



4.3 Typical entanglement spectrum 



Once the exact solution is found one can study any quantity of interest. Using general results of elementary algebra, one 
can extract plenty of information from the coefficients of the Laguerre polynomial. For example, in order to compute 

the purity, Eq. (5|, note that ^i) = J2i + '^J2i<j ^^^ji t^^^t is 



T^AB 



{N-l){M-2) N + M-2 



N{M - 1) 



N{M - 1) 



(41) 



which is exactly the same as Eq. (31 1, to be compared with the average value computed by Lubkin, Eq. (14). For large 
N ^ TTAB — 0{1/M) while the difference between mean and mode is 0{1/NM). Similarly, one can readily compute the 
elementary symmetric invariants 



E 



.A, 



jl<j2<---<jk 

m{M -1)1 



k\{N - k)\(M - k - l)\ 
As an example, the typical determinant of the reduced state pA is 

{M-l)\ 



CN-k 
CN 



[NiM~l)Y 



sn 



(A) = J|A, = det (pa) 



(M - N -1) 



[N{M ~ 1)] 



-N 



(42) 



(43) 



Let us conclude this section with a remark. A given property is said to be "typical" if it holds with overwhelming 
probability. In the previous section we just found the most probable spectrum of the reduced state pA- The next step 
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whould be to show that the probabihty concentrates around the most probable value. Some authors [221123] achieved 
some exact results on the concentration of the spectrum of the reduced density matrix when both subsystem are large. 
We will give here a heuristic justification of this assertion based on a "second moment bound" . In what follow we will 
specialize our discussion to highly unbalanced bipartitions and a large-size limit. Let us suppose that subsystem A is 
very small compared with B so that we can think at B as an environment. One expects that when the dimension of 
the environment M is large, the eigenvalues fluctuate around 1/N. Thus we make the educated guess 



N Af" 



(44) 



The value of the variance of the random variable in the unbiased case, Eq. (12), suggests that a = 1/2, and the 
saddle point equations in terms of the e^'s read, for large values of M 



2E 



1 



+ nVm - N^Ei- ^^/aI ^0 , yie{i,...,N} 



0. 



(45) 



We can use Stieltjes's trick to handle the above simplified saddle point equations. In this framework, the first N 
equations in ( 45 ) are equivalent to the single differential equation: 



g"{x) - {N^x - nVm + ^VM)g'{x) + N^gix) = 



(46) 



whose polynomial solution Hj^ is the Hermite polynomial of degree N [2T]. Therefore, the most probable eigenvalues 
satisfy 



H 



N 



V2 







(47) 



This approximated solution is useful to find a second-moment bound with a Gaussian approximation. The most 
probable value is typical (and close to the average) if the pdf is sharply peaked. The width of the peak is encoded in 
the Hessian H{\) of the energy, Eq. (23), evaluated at the minima. In the Gaussian approximation, the proper values 
/li's of H give the widths (cr~ = hi) of the pdf around its maximum. In order to extract the correct scaling, we can 
consider just the trace of the Hessian and use the solution of the simplified saddle point equation, Eq. (47), to obtain 



trff < N^{M ~N) + 2N{N - l)M 



(48) 



Since ii H = 1/cf < N/ ■ , we readily find an upper bound on the width of the Gaussian. 



4.4 Thermodynamic limit 

In this section we want to use our solution to investigate the large size limit of two macroscopic partitions. As in 
section [2] it is convenient to introduce the parameter [i such that ^iN = M — N . When N,M ^ +00 with jj, fixed and 
finite, the most probable purity reads 

M + N-2 N{2 + n}-2 l{2 + ^l) 

a value obtained in [B] by more sophisticated methods. 

Similarly, one can compute the quantities trp^, and then derive the Renyi's entropies. A very inexpensive method 



is to write ^^X^ in terms of elementary symmetric polynomials, and then use Eq. (42). The above procedure can 
be easily implemented by the most common symbolic manipulating softwares {Mathematica provides the suitable 
function SymmetricReduction [] for expanding any symmetric polynomial). We give a list of the first five traces in 
the thermodinamic limit, M — > +00 with /i fixed: 



, _ {2 + ^l) 1^/1 



^^•^^= (i+> ' 
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4 14 + 21^ + 9^l^ + /i3 1 / 1 ^ 

ti-PA^ ^^1^ TH+OKtI ' (52) 



(1 + ^)3 iV3 V^"* 

trPA- , . X4 T?4+0 h^ . (53) 



(l + /i)4 iV4 V^'^ 

The same can be done for the iV-degree elementary invariant, that is detp^- For example, for a balanced bipartition 
we have 

ATI 

detPA-^. (54) 



5 Typical states of fixed entanglement 



So far, we have dealt with unbiased pure states. We now look for the most probable eigenvalues sampled on isopurity 
manifolds, i.e. the most probable value of the pdf fN.Aii^), given in ([9]), on the manifolds 



(55) 



which geometrically are given by the intersection of the N — 1-dimensional sphere of radius ^/ttab with the simplex 
Z\jv-i- The new constraint, — '^ab, enables us to compute immediately the second moment of Ai, that is, the 

width of the pdf restricted on the isopurity manifolds. From the permutation invariance of fN,M- 



and then, the variance is: 



1^' 



l<i< N 



{{^')) = (((A.) - A.)^) = (A,:^) - (A.)^ = - 1), l<^<N 



(56) 



(57) 



where, in the last equality the constraint on the simplex was used. 

By introducing a second Lagrange multiplier rj to take into account the new constraint, one has to find the minima 
of the {N + 2)-variable function 



^tot (A, t log I A. - A, I - {M -N)Y, log A, - ^(1 - ^ " ^(""ab - E ^^)- 



(58) 



The saddle point equations ( 24 ) modify into 



-277A, + 2^ 



A,; — Ao 



M -N 

A,; 



1< i < TV, 



(59) 



As done before, by multiplying the first N equations of system ( 59 ) by A^ and taking the sum we find 

^ = N{M - 1) - 2t]7tab ■ 



(60) 



As pointed out before, to find the most probable spectrum is analogous to the problem of finding the equilibrium 
positions of interacting charges in an external field. The charges repel electrostatically via a 2D-Coulomb potential 
in an external potential that now includes a new term due to rj ^ 



fc^Lt {x) — ~a\n X + £_x + rjx"^ 



(61) 



where a = M — N > Q. Kt the origin there is a charge a that repels the other charges. Moreover, the unit charges 
experience a constant force whose direction is opposite to the sign of ^, and the Lagrange multiplier rj plays the role 
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of the elastic constant of a harmonic potential tjx^. Thus, we have reduced the problem to that of finding, among all 
configurations of the Coulomb gas with constrained center of mass and momentum of inertia, the one with smallest 
electrostatic energy. Another, equivalent, point of view, is to look at rj as an inverse temperature which fixes the energy 
ttab of the system, and will be discussed in section [6| 

As discussed in the previous section, when 77 = the Lagrange multiplier ^ has to be positive in order to have an 
energy minimum in (0, 1). The minimum is at Xc = Q^/C- Since we expect that this critical point is close to 1/N we 
can immediately guess ^ ^ NM in the large M limit. When rj the scenario becomes more interesting. Again the 
N charges will arrange themselves in the external potential in a configuration of minimum energy. In the following 
analysis, we will consider the balanced situation a ~ M — N = 0. 



5.1 Balanced bipartition 



The balanced case, AI — N, is exactly solvable in the context of orthogonal polynomials. The saddle point equa- 
tions (|59|) specialize to 

- ^ = , l<i<N, 



-2t]X, -I- 


2E 












T^AB 



Aj - A, 



(62) 



By summing over i and by using Eq. ( 60 ) with M = N we get 



2r] 

'TV 



T^AB 



1 

N 



N{N - 1) 
2^ 



(63) 



The Stieltjes method leads us to the differential equation 

g"{x) - (277X + Ogi^) + 2T]Ng{x) 



(64) 



4 



rjx 1 . We fix the Lagrange multiplier ^ 



whose polynomial solution is the Hermite polynomial [21j Hn 

—2ri/N in order to satisfy the unit-trace condition. Then, the most probable eigenvalues are the N solutions of the 
following polynomial equation: 

1 

X 

N , 



= 



(65) 



Moreover, from Eq. (63), we find how 77 labels different isopurity manifolds 

N^{N -1) 



77 = 



2iNnAB - 1) 



(66) 



Therefore, since 1/iV < ttab < 1, one gets that rj > N'^/2, the minimum being attained for separable states, tt^b — Ij 
while 77 — +00 for ttab — > 1 /N. 



However, if 77 is too small the solution to (65) ceases to be physical, since one or more eigenvalues become negative. 
That means that the most probable spectrum belongs to the boundary of the simplex, and it is no longer a critical 
point of the total energy. Thus, for all TV, there exists a threshold value 77+ > 7V^/2 below which the solution suddenly 
ceases to exist. Now, recall that, denoting by ctat the smallest zero of the Hermite polynomial Hn{x), one gets \2T\ 



lim 



=-V2 
'N 



(67) 



For large positive values of i] we expect to obtain very mixed states with ttab = 0{1/N). Then, in the large N limit, 
from Eq. ( 66 ) the correct scaling for the Lagrange multiplier 77 is 



77 = /37V^ /3 > 0. 



(68) 



Using the above asympotic property of the largest zero of the Hermite polynomials, we readily find Nxj^ — > 1 — y^^T/^i 
whence the critical value is /3+ —2, and the corresponding critical values of the purity tt^^ = 5/4 A^, in agreement 
with refs. [6l[7]. The typical eigenvalues of the fixed-purity ensemble are the zeros of Eq. (65 1, whenever ttab < 5/4A. 
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6 Canonical Ensemble and Partition Function 

The above approach is based on the microcanonical ensemble, in which the purity of distinct manifolds is fixed. A 
different way to proceed is to fix the average purity by introducing a partition function and reformulate the problem in 
terms of a classical canonical ensemble. The main quantity we are interested in is the local purity ttab as a measure of 
the bipartite entanglement between balanced bipartitions. This quantity will play the role of energy in the statistical 
mechanical approach. 

Let us clarify the rationale behind our analysis. Although our interest is focused on the microcanonical features 
of the system, namely on "isoentangled" manifolds, we find it convenient to define a canonical ensemble and a tem- 
perature. This makes the analysis easier to handle and it is based on the equivalence -largely used in the statistical 
mechanical description of large systems- between the microcanonical ensemble (in which energy is fixed) and the 
canonical ensemble (in which temperature is fixed). 

The inverse temperature /3 is a Lagrange multiplier for the optimization problem. It is the variable that is naturally 
conjugate to ttab- P fixes, with an uncertainty that becomes smaller for a larger system, the level of the purity of 
the subset of vectors under consideration, and thus an isoentangled manifold. The use of a temperature is a common 
expedient in minimization problems that can be recast in terms of classical statistical mechanics. We notice that this 
approach has been fruitful in several other context such as the analysis of mixed states |24j and multipartite pure 
entangled states P51E51[?7] where entanglement exhibits the phenomenon of frustration [55] . 

Therefore, in order to study the typical properties of a large bipartite quantum system, we introduce a partition 
function from which all thermodynamic quantities can be computed 

Zab = f exp {-I3N^7:ab) dv{\) , (69) 



Af-l 



where dv{X) = f]^^pf{\)d\, with pdf Q. Here, N^ttab plays the role of "energy" while ^ is a "temperature" that 
control energy, that is entanglement. The factor is chosen in order to make the energy an extensive quantity, 
N^TTAB = O (N^), Since {ttab) = O (l/N). 

For large N we look at the maximum of the integrand, that is the maximum of 



Zab = / exp(-7VV(A,C,/3))dA, 



(70) 



where 

ViX, C, /3) = /3A^ ^ - A ^ In I A, - Afe I + C( E ^fc - 1) ■ (^1) 

j j<k k 



As before, the Lagrange multiplier Q fixes the normalization constraint and allows one to extend the integration from 

il = f3N\ ^ = QN\ (72) 



the simplex Am-i to all positive values M?. Notice that, by setting 



in agreement with Eqs. (68) and (63 1, we get 

iW(A, C, P) = Etot{\ v) + mAB, (73) 



where i^tot is the energy (581 for a balanced bipartition, N = M, thus establishing the equivalence between the two 
approaches of Sec. |6] and Sec. 5.1 Indeed, in the thermodynamic limit, the main contribution to the integral of the 



partition function is given by the maximum of its integrand. The standard way to solve the problem is to apply a 
saddle point method, i.e. to look for the stationary point of V, that is the stationary point of i?tot, which gives again 



Eqs. (62). Notice, however, that in the canonical framework, the last equation of ( |62[ ) is not viewed as a constraint, 
but rather as a relation between temperature and average purity. 



Moreover, in a thermodynamic analysis it could be convenient to rewrite the saddle point equations (62 1 in the 
continuous limit. To this end, by recalling that all the eigenvalues A^ are of order 0{1/N), we introduce the empirical 
distribution 

j 

that in the limit of large N can be approximated by a continuous density function. Obviously, we get 

cr(A)dA = l, (75) 
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N = M 



Max Entangled State 
^AB = 1/N 

( /3 = +CO ) 



^AB = 5/4N 



^AB = 1 / (N-1) 



Phase Transition 

( /3 = /!+ ) 



A 



Fig. 4. Left: Density of the eigenvalues for /? = 2, 4, and 10: analytic solution, Eq. (78 1 and numerical results (zeros of Eq. 
(65 1). In the temperature range j3 £ [2, +cx3) the solution is given by the semicircle law. Right: Density plot of the pdf fN,N{X), 
Eq. in the simplex of eigenvalues of the density matrices obtained by par tial tracing in the balanced case A'' = M = 3. The 



green line shows the locus of the most probable spectra (solution of Eq. (65l) as /3 = rj/N varies. For j3 — >■ +oo the spectrum 



tends toward the maximally mixed one {1/N, . . . , 1/N), corresponding to a maximally entangled pure state of the total system. 
As P decreases, the typical state becomes less entangled (reduced density matrix less mixed) and reaches the boundary of the 
simplex for j3 — j3+. In the thermodynamic limit, at /3+ = 2 a phase transition occurs, one eigenvalue vanishes, and a new 
equilibrium configuration (Marcenko-Pastur distribution) takes place instead of Wigner's semicircle law. 



and by making use of ( 74 1 , the saddle point equations (|62| can be easily rewritten as 

J X- fi 2 



A(t(A) dX = 1, 

J A^ct(A) dX = TTAB , 
where /x = A^A^, / is the Cauchy principal value integral, and 

TTAB = NtTAB- 

The first equation is a singular Fredholm equation of the first kind, known as Tricomi's equation [5S 
lies on a compact interval (depending on and can be given explicitly 29_. One obtains for /? > 
semicircle law 

a(A) = ^V(A-A_)(A+-A), 

TT 



(76) 



where (see ref. [7]) 



A+ 1± 




1 

2P 



(77) 

Its solution 
= 2 Wigner's 

(78) 



(79) 



In Fig. |4]we plot a{X) for several values of /3. Observe that as /3 becomes larger the distribution becomes increasingly 
peaked around 1. This means that in (74| all the eigenvalues tend to 1/iV: for temperatures 1//3 close to zero the 
quantum state becomes maximally entangled. 

At higher temperatures < /3 < the solution acquires a different physiognomy, namely the Marcenko-Pastur 
law [3D]. The change from semicircle to Marcenko-Pastur, a breaking of the Z2 symmetry of the solution about 1, is 
known as the Wigner-to-Wishart second order phase transition [5]. We just mention that, for positive /3's, lowering the 
temperature will make the state nearest to the separable ones. As a consequence, among all the isopurity manifolds 
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the one corresponding to /? = 0, namely the unbiased ensemble, exhibits the highest purity: 



a(A) 




4- A 
A 



T^AB 



^ 2 . 



(80) 



In order to obtain more separable typical states, the sampling has to be done with negative /3's. By lowering /3 
below zero, in the framework of the partition function one finds that cr(A) experiences many others phase transitions 
[7] . The analytic continuation of the stable solution for positive temperatures emanates a metastable branch in which 
all eigenvalues remain 0{1/N) (and so does the purity) even though the temperature can be negative. By constrast, 
a new stable solution emerges by allowing the largest eigenvalue to evaporate, leaving the Marcenko-Pastur sea of the 
other eigenvalues 0{1/N) and becoming of order of unity. This is a signature of separability: negative temperatures 
explores regions of the pure state space mostly populated by separable states. 



7 Conclusions 

In this paper we have analyzed the properties of the typical spectrum of the reduced density matrices obtained from a 
suitable ensemble of pure states. We have found that the most probable eigenvalues of the density matrix coincide with 
the equilibrium positions of charges on a line with interaction forces arising from a logarithmic potential. In particular, 
for an unbiased ensemble, the solution of the problem can be obtained from the zeroes of Laguerre polynomials. In 
the case of a fixed value of purity, the ensemble is no longer unbiased and the charges of the model are affected by 
a supplementary quadratic term in the external potential. This leads to a modified solution and the most probable 
eigenvalues are given by the zeroes of Hermite polynomials. Finally, in the case of balanced bipartitions, we have used 
an approach based on a temperature and a partition function. In this case we have obtained Wigner's semicircle law, 
and a second order phase transition to a Marcenko-Pastur law. 
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